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OF  TRANSCENDENTAL  FUNCTIONS 


by 

Jerry  L.  Fields 
and 

Yudell  L.  Luke 


Summary  and  Introduction 


In  a previous  report  [l]i/,  j^t  was  shown  that  a large  class  of 
transcendental  functions  could  be  represented  by  a rational  function,  the 
ratio  of  two  polynomials,  together  with  a remainder  term.  a More  specifically, 
if  E(z)  is  defined  by  a Laplace  Integral,  then 


JS(t)  - ^/fn  +R„  , Rn  * Vfn 


where  Cpn  and  Fn  are  polynomials  of  degree  n and/Rn  is  the  remainder. 

In  Section  I it  is  shown  that  this  can  be  considered^  regular  summation 
technique . 


By  proper  choice  of  some  parameters,  fn  is  a polynomial  of 

hypergeometric  form.  If,  except  for  a multiplicative  factor,  E(s)  has  a 
generalized  hypergeo me trie  series  representation,  then  Fn  is  a sum  of  n 
generalized  fcypergeoraetric  series . 

^ Only  in  certain  cases  could  it  be  shown  that  the  rational  represen- 
tation converges,  i.s. , for-  a — fflmeA, — Us — - 0 One  of  the  prin- 

n-*oo  ^ 

cipal  difficulties  lies  in  the  fact  that  asymptotic  expansions  of  a certain 
class  of  hypergeometric  polynomials  were  known  only  for  a few  special  cases, 
the  classical  Jacobi  polynomials,  for  example.  Asymptotic  expansions  of  this 
class  of  hypergeometrlc  polynomials  now  exist,  eee  [t^j>  and  the  convergence 
question  can  be  further  explored. 


\ 


I 7 Numbers  in  square  brackets  pertain  to  references  at  end  of  report. 


- 1 


The  structure  of  the  numerator  of  the  error  term,  i.e.,  Pn  is 

complicated  and  only  in  a few  special  cases  has  it  been  possible  to  represent 
it  in  a useful  fora.  This  dravback  is  partly  alleviated  by  the  fact  that  one 
of  the  above- mentioned  special  cases  includes  the  Whittaker  functions  and 
their  natural  generalisations.  For  those  cases  where  FQ  can  be  put  into  a 
useful  form,  we  prove  convergence  of  the  rational  representation. 

In  the  case  of  particular  interest,  i.e.,  the  Whittaker  functions, 

we  represent  Fn  as  an  integral  and  then  find  a suitable  bound  for  it.  This 

coupled  with  an  asymptotic  estimate  of  fQ  enables  us  to  prove  that 

lim  Rn  « 0 and  also  to  obtain  an  asysptotic  bound  for  the  error.  This  is 
n — » o© 

the  essence  of  Section  I.  In  Section  II,  numerics  are  presented  to  display 
the  effectiveness  of  this  bound. 


1.  Convergence  of  the  Rational  Approximations 


In  [l],  it  was  shown  how  to  formally  generate  a rather  general 
sequence  of  rational  approximations  for  the  generalised  hyper geometric 
function  [j] 


fi(t) 


p+lpq 


• • . , 

. . •* 


(1.1) 


where 


«FnW  * ofn  (*1*"’*  °‘» 


f t* ( 

k-°  ¥ ( ^)k 

t»l 


kt 


(1.2) 


Throughout  this  paper  we  employ  a contracted  notation  and  write 


2 


mPn<*>  - mV 


«i 

ft* 


oo 


(1.3) 


Thus  ( <xm)  is  interpreted  as  71  ( c*t)  , and  a similar  remark  holds  for 
k t«l  k 

( /°n).  • We  assume  that  none  of  the  o^' s , or  /°t'B  are  zero  or  a 

negative  integer.  Also,  it  is  assumed  that  the  difference  of  any  numerator 
parameter  , and  any  denominator  parameter  , is  never  equal  to 

zero.  Empty  terms  in  expression  such  as  (1.3)  are  replaced  by  unity.  Two 
distinct  sequences  of  rational  approximations  were  developed  in  the  above 
paper  [l] . They  are  as  follows. 


Type  I or  Homogeneous  Case 


EU)  - ♦ R„U,*)  , RnU,/>  - * • <!•*> 


r „[r](/> - Th-e/  • f0[°]  ' ) 


(1.5) 


(1.6) 
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i ' 


oo 


'>  ■ TiTT g S ■%} 


^^r+k+l^*^ 


r+1 


r+k+1 


(r+k+1 )l 


dt 


(1.7) 


oo 


./e-ztt^-i^an_ /XJ.AV  (3fpWi-i(q'+ra)k(Xt)r+1at 

0 ic^O  \*  (>•_)  . . (r+k+1 )'. 


T(cr) 


^ ^<1  r+k+1 


(1.8) 


^ an-k(°rir+l)ic(  orp+r+l)  k 
/ ^0  (2-tr)^  />q+r+l) 


* ft  & teM  /°q+1,r 


k/X^T  \ 

v U / * 

(1.9) 


p q 

idiere  rx^  = IT  oc^.  > and  = TT  ^ • Here  the  coefficients  a\»  are  as 
t*l  t=l 

yet  arbitrary,  and  * is  unrestricted  except  that  | i/z\  <1  . 


Type  II  or  Mon-Homogeneous 


B(z)  = vV  + Rn(z»  * > * Rn(z>  * ) = F-n(%2/  * ] 

rn'  • f fn(  i ) 


(1.10) 


(1.11) 


M 


fnLrJ(*)  ■ ^ »n-k  J*  - fJ^OT)  - f or ) 


(1.12) 
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( /°n)^*v(r+k) 


' s are  unspecified  and  V is  subject  only  to  the 


The  difference  between  the  two  sequences  is  that  when  z = t,  in 
the  Type  I representation,  pn(0)  ^ 0 , whereas  for  the  Type  II  representa- 
tion under  the  same  condition,  9?n(0)  = 0 . For  seme  particular  case  one 
of  the  representations  may  be  more  desirable  than  the  other.  Also,  it  some- 
times happens  that  one  sequence  approaches  the  desired  limit  monotonically 
from  above,  and  the  other  monotonically  from  below,  thus  yielding  rational 
inequalities  for  the  functions  in  question.  The  above  formal  sequences  were 
derived,  strictly  speaking,  with  the  restrictions  p<q  and  Re(z)>0  if 
p<q  , while  Re(z)>Re(x)  if  p=q  . However,  these  restrictions  may  be 
weakened  whenever  the  resulting  expressions  are  defined  and  make  sense. 


The  merit  of  [lj  lie/s  in  the  fact  that  it  motivates  how  the  arbitrary 
coefficients  should  be  chosen  and  gives  a closed  form  expression  for 

the  error.  The  a^'s  were  chosen  in  [l]  to  agree  under  certain  circum- 
stances,at  least,  with  the  rational  approximations  given  by  application  of 
Lanczos*  r-method  (see  [4])  to  the  linear  differential  equation  satisfied 
by  E(z)  (see  [3]).  The  choices  made  for  the  a*  were 


^-k 


■W 


if 


r b l 


(1.17) 


(<Vl)k(<r+Dk  xk 


* 


where 


cnk 


(-n)k(n+<y+(3-n)k 

((5+l)k  ki 


Type  II  - Non-Homogeneous 


Cnk( 


(1.18) 


and  is  defined  as  above.  With  this  selection  of  the  coefficients  ak  , 

we  show  that  the  rational  sequences  defined  by  (1.4)  - (1.9)  and  (1.10)  - 
(1.15)  converge  for  two  general  classes  of  functions  defined  as  follows: 

Class  I is  composed  of  those  hypergeometric  functions  such  that  piq  in 
(l.l)  and  Class  II  is  made  up  of  the  hypergeometric  functions  such  that 
p * q+1  in  (l.l).  Class  II  Includes  the  Whittaker  functions,  modified 
Bessel  functions,  the  Weber  Parabolic  Cylinder  functions,  and  various  special 
cases  thereof.  For  completeness,  we  define  a Class  III  composed  of  all 
hypergeometric  functions  such  that  pi  q+1  in  (l.l). 

Considering  E(z)  as  a formal  infinite  series,  we  define  the 
partial  sum  of  this  series  by 


_ ^ (o-)r(«p)r  (A/z)i 


(1.19) 


Rearrangement  of  (1.5)  and  (l.ll)  gives 


- T H-njrk  *k<*> 

K=a 


- 6 


(1.20) 


Mkwuiiaii 


where  a * 0 for  the  Type  I or  Homogeneous  Case, 

a a 1 for  the  Type  II  or  Non-Homogeneous  Case. 

Thus  our  rational  approximations  may  be  considered  as  a summation  technique. 

A method  of  suomability  is  said  to  be  regular  if  it  sums  a convergent  series 
to  its  ordinary  sum.  Vte  now  show  that  our  economization  process  is  a regular 
summability  process. 


Theorem.  If 


(1) 

* ) = an-k  ^ 

= It-  an-k  ^ 

k*a 

k=0 

(2) 

Lim 

P,,(z)  = E(z)  , for  fixed  z , 

k— ► oo 

(3) 

an.k>0  , / >0  , and 

(4) 

Lim 

T1—+00 

tnU  ) = °°  * fa+1(  * ): 

>'n<*> 

for  fixed  if  , 

then 

Lim  |pnU  / )/fn(  { ) 

n— »oo|  u 

- E(z) 

= 0 for  fixed  ^ and 

Proof. 


Given  € >0  , there  exists  a positive  integer  N (usually  dependent 
on  z)  such  that  n>N  implies  |Pk(z)  - E(z)|  < e/3  . Then 


i )/fn(  i ) 


[Pk(z)  - E(z)  - 


an-k  >k  [^(z)  " E(z)  ] - a an  [c/3  + E(z)] 

7JT) 


- 7 - 


Thus,  for  n >H  , 


<pn(*,  jr  >/fn(  jr)-E<Z) 


ian-k  <k|Pk(*)-*(0-£/s|  j - a «a(e/3.|«.)|  ) 


+ €/3 


fn<») 


Since  N is  fixed,  n can  be  chosen  large  enough  such  that  the  right  hand 
side  of  the  above  Is  less  than  €"  . This  proves  the  theorem.  From  the 
results  of  [2]  , and  (1.16)  - (1.16),  we  see  that  convergence  for  Class  I of 
the  hypergeometric  functions  is  proved,  when  A 20  . 

Ms  now  consider  Class  II,  i.e.,  the  situation  where  p = q+1  . 
Treating  the  non-homogeneous  or  Type  II  < ase  only,  for  p = 1 and  q = 0 , 
from  (1.15)  and  (1.16)  we  get 


r(a-)  0 


00 

yVztt<r-l  21  W jr/z)k  3F2  A*i+k,<r+k,l 
3 k=0  \ k+l,cr 


k+l,o- 


At  dt 


(1.21) 


By  use  of  the  relationships,  see  [3] 


p+lIVl  lz)  * r£$7^-2)  / (£l“)  du 


Re  0 , Re 


(1.22) 


and 


(1.23) 


we  can  write 


• 8 • 


oo  1 


r“(8’ * 1 ' r(<r-i)r(i-«1)r(«1)  / / / 


/ a-ZttO~  -lu^l-l(l.u)- ^ly^"2 

(1-t  X u) 


x (&£))—•>  - 

Re  X< 0;  Re(z)>0,  0<Re«i<l,  Re<r>l 


(1.24) 


Defining 


Mil  « max 
71  0^  vil 


3P2  ( ~n>MOf+P  +1*1 

\ P+1,  <*i 


(1.25) 


and  noticing  that 

1 1-t  Xu  I ^ 1,  Oitico,  0£u£l  , 


0 £ 


“(-rK&rf) 


(1.26) 


it  is  easily  seen  by  direct  confutation  that 


Pn(«#  jrJiMf1 . 


(1.27) 


By  repeated  use  of  (1.22)  and  the  same  method  of  proof  as  above,  it  can  be 
shown  that 


11  _ max 

oi  vi 


MiA 


(1.28) 
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Prom  [2),  if  (2t»>+a  + 0+3/2)  >0  , where 


2<*> 


i/2  + ^ <*  - (l  +§♦  ^ Oti^) 

i«i  V i«i  y 


(1.29) 


*$L 


I-1  (/® -d.)P  (2«>+ff+  j5+3/2)  ' J 


N2  ■ n(n+o<+ £+1) 


(1.30) 


Thus  Fn(z,  /)  is  bounded  by  e term  of  algebraic  order  in  n . 
(*1  the  other  hand,  by  (1.16),  (1.12)  and  a result  given  in  [2], 


- 10  - 


where 


N2  ■ n(n+<x+(S+l) 

do  - 2«»+l/a-<r 

3 

X - - X , Re  (X)iO 


<P(-y/x)  - + - (oc+0+2/3+2t3)( y/x ) 

(Dl-I^)  , 

+ ■ ■ (2Di.+I^-l)  +E2-Ei-2/9 

w 

q o+3 

°L  * 1 + Z (^V1)  * d2  ■ Z <*1  * <Ve  ■ * wq+3  - 1+  P 


*0.  - Dl-1  + £ f;  (^n-D(/,j-l) 

n-2  J-l 

q+3  n-1 

% - Z Z 

n-2  J-l 


(1.32) 


Since  f n( % ) is  of  exponential  order  in  n , 


Fn(sj  2T  ) 


(1.33) 


Similarly  in  the  homogeneous  or  Type  I case,  one  can  show 


where 


max 

Oivil 


<r+  l 


max 

0<.  w<  1 


n.,P_.? / -n,  n+  or  + 0 +1,  2 I w\ 


CT*  - 1 


Again,  by  [2]  , it  (2n+<X+  0+3/2)  >0  , where 


2n 


i + qr  0 + £ A - 2l  oc( 


~ «o+l 
1-1  i-1 


/I 


§-*- 


♦ 1 fl 

i«l 


i * 


<1+1 

H aq+l 

1*1 


<r«  1 


\h 


r(l+/C5l)r(2n+or  + ^+3/2)  1 Wj 


/I 


P(  /Oq)r (2/1+^+ (3+3/2)  l \N/J 


cr  ■ l 


> . 


(1.35) 


(1.36) 


(1.37) 
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M2  ■ n(n+cf+p+l) 


Thus,  by  (1*34)  and  (1.37),  Fn(z,  $)  is  again  bounded  by  a term  of  algebraic 
order  in  n . Combining  thi3  with  (l.ol),  we  again  hove 


I 


11m 
n— * 


ooRnU>*) 


lim  Fn(z,  *)  _ 

_ k ■■I—."  ■ m 0 

n~*°°  fn(y) 


(1.36) 


It  calls  for  remark  that  the  bound  for  the  integral  in  (1.24)  was 
obtained  by  ignoring  the  oscillatory  nature  of  the  hypergeometrlc  function 
in  its  integrand  over  the  entire  region  of  integration.  Thus,  the  asymptotic 
bound  for  the  remainder  is  conservative  as  the  rational  approximations  con- 
verge much  more  rapidly  than  indicated  by  (1.26).  Realistic  estimates  of 
Fn(z,  ti ) Qn<1  80  of  Rn(z>  f ) seem  to  be  much  more  difficult.  However, 
though  (1.26)  is  quite  rough,  it  is  easy  to  use  and  not  excessively 
misleading. 


2.  Numerics 


To  give  a qualitative  idea  of  the  effectiveness  of  the  bound 

Mq+{  for  F0(z,  jf ) , we  consider  an  example  given  in  [l]  (Table  (6.2)).  The 
example  is  connected  with  the  modified  Bessel  function  Ky  (z)  . 


Kv  (z)  - e"*(r/2 z)^2S(z) 

E(z)  - 2F0(l/2  -V,  1/2  +V|  - l/2z) 

Using  the  non-homogeneous  or  Type  II  approach, 


(2.1) 


13  - 


rn(  ) - 3f3  (ma> n+  or+  Q +1*1  I -2  v) 
n 3 5 \^+l,l/2’V,l/2+y  I ) 


U*>*)  - (1/4- v2)  Z 


Cn  yTy/z)1 


r-0  (l/2+V+r)(l/2-V+r) 

\ 

XP  ( -n+r,n+<*+ 0+l+r,l  I 

3 3 @ +l+r,3/2+z>  +r,3/2-  V +r  | 


(2.2) 


- eq(z>  jr ) + Rn(z>  ir) 


If  oCs  ^ ■ -l/2  , v = 0 , X m z > and  denoting  asymptotic  values  as 
computed  from  (1.30)  and  (1.31)  by  , we  have  to  four,  significant  figures, 


f3(z) 


'f*(z) 


*3(2,  jQ 

0.1670 

0.3977 

0.3605 

0.0— 

-0.1843 

-0.4053 


Rj(^  if) 


0.5 

298.3 

295.5 

5.598 

X 

IO-4 

1.0 

1 

265. 

1 

276. 

3.145 

X 

10-4 

2.0 

6 

515. 

7 

007. 

5.84 

X 

10-5 

5.0 

72 

400. 

121 

000. 

0.  — 

10-6 

7.0 

184 

300. 

507 

800. 

—1. 

X 

10.0 

506 

600. 

3 994 

000. 

-8. 

X 

10-7 

/~f3(z) 


3.188  x 10“2 
7.382  x 10*3 
1.344  x 10-3 
7.785  x 10-5 
1.855  x 10“5 
2.359  x 10“6 


(2.3) 
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f4(g)  R4(t,  t) 


0.5 

1 206. 

1 198. 

1.170  x 10*4 

1.0 

7 169. 

7 217. 

-6.1 

x 10-6 

2.0 

55  840. 

61  170. 

-5.4 

x 10*6 

5.0 

1 208  000. 

1 754  000. 

-2. 

x 10-7 

7.0 

4 052  000. 

8 491  000. 

—1. 

x 10-7 

10.0 

15  140  000. 

68  760  000. 

-1. 

x 10-7 

z 

F4(*,n 

Mi11 

II 

'''Mj11/ r^f4(s) 

0.5 

0.1411 

12.58  12.57 

1.049  x 10*2 

1.0 

-0.04375 

1.742  x 10*5 

2.0 

-0.5015 

2.055  x 10-4 

5.0 

-0.2417 

7.166  x lO*6 

7.0 

•0.4052 

i . 

1.480  x 10"6 

10.0 

-1.514 

\ 

f s 

f 

1.828  x 10-7 

(£.4) 


The  above  example  shows  that  our  error  analysis  is  quite  conserva- 
tive, as  expected.  For  the  present,  a pragmatic  view  should  be  taken  con- 
cerning the  error.  That  is,  IT  more  precise  information  is  required,  in 
those  cases  where  convergence  is  assured,  one  should  compute  for  successive 
values  of  n and  accept  the  conmon  digits  as  correct.  To  illustrate,  in  the 
above  example  for  K^s)  , if 

Ej(0.5)  • 0.85952 

E4(0.5)  - 0.85976 

and  it  is  quite  reasonable  to  say 

E(0.5)  - 0.860 


with  a possible  error  of  one  unit  in  the  third  decimal  place 
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